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Although the Lepton Flavor Violating (LFV) decay fi — > e+7 is forbidden in the Standard 
Model (SM), it can take place within various theories beyond the SM. If the branching ratio of 
this decay saturates its present bound [i.e., Br(pi+ — > e+7) ~ 10 -11 ], the forthcoming experiments 
can measure the branching ratio with high precision and consequently yield information on the 
sources of LFV. In this letter, we show that for polarized fi + , by studying the angular distribution 
of the transversely polarized positron and linearly polarized photon we can derive information on 
the CP-violating sources beyond those in the SM. We also study the angular distribution of the 
final particles in the decay fx + — » e+e~e+ where e+ is defined to be the more energetic positron. 
We show that transversely polarized e+ can provide information on a certain combination of the 
■ CP- violating phases of the underlying theory which would be lost by averaging over the spin of e+. 

in 

PACS numbers: 11.30.Hv, 13.35.Bv 
^■/■^ ■ Keywords: Lepton Flavor, Muon decay, CP-violation, angular distribution 

> 

O 



o 
o 



o 
o 

Oh! 



X 



I. INTRODUCTION 



While the Standard Model (SM) preserves the lepton flavor, its various extensions such as supersymmetry or large 
extra dimensions can lead to Lepton Flavor Violating (LFV) rare decays /j+ — * e + 7 and /i+ — + e+e~e+ detectable in 
the forthcoming experiments. The present experimental bounds on the branching ratios of these processes are [l[ 



Bt(p+ -> e+ 7 ) < 1.2 x 1CT 11 Br(^+ -> e+e+e") < 1.0 x 10" 12 at 90% C.L. 

The MEG experiment at PSI 0, which is under construction, will be able to probe Br(/i+ — > e+7) down to 10 . 
Thus, if this branching ratio saturates the present bound (i.e., Br(/i+ — > e+7) ~ 10 -11 ) the future searches will 
enjoy high statistics and can make precise measurement limited only by systematics. Moreover, since the muons 
are produced by decay of stopped pions (at rest), they will be almost 100% polarized. Thus, studying the angular 
distribution of the final positrons, we can learn about phenomena such as parity violation, through which more 
information on the sources of LFV can be extracted [H . Among the various extensions of the SM that can give rise to 
lepton flavor violating phenomena, in the literature the Minimal Supersymmetric Standard Model (MSSM) and large 
extra dimensions have received particular attention. It is well-known that in both cases integrating out the heavy 
states of the model, the LFV Lagrangian responsible for fi — ► e7 can be written as 

L = A R fi R cr^e L F^ + A L fl L a^e R F^ + A R e L a^fi R F^ + A* L e R a^fi L F^ (1) 

where — 5 [7**, 7"] and F^ is the photon field strength: F M „ = d^e^ — d u e^. Although Lagrangian in Eq. JTJ is 
not the most general form of the effective LFV Lagrangian, throughout this paper we consider only these terms. In 
the appendix, we consider a more general form of the effective Lagrangian and show that (|TJ) is indeed the dominant 
part. 

It can be shown that averaging over the spins of the final particles, Lagrangian © yields H 

dV Kc2e +l) = [\ A R\ 2 {l-V a cosQ) + \A L \ 2 (l + P u cos6)} , (2) 
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where is the angle between the momentum of the positron and the spin of the muon and P M is the polarization 
of the muon. Notice that integrating over cos0, we arrive at T(p + — > e +r y) = (rrv^/ Att)(\Al\ 2 + |A^| 2 ). Thus, by 
measuring the total decay rate to e + 7, we can only measure \ Al\ 2 + |^4i?.| 2 - However, Eq. ([2]) shows that by studying 
the angular distribution of the final particles with a moderate angular resolution, |^4/j| 2 and \Al\ 2 can be separately 
derived. Information on |Ai,| 2 /|Ar| 2 can be translated into information on the sources of lepton flavor violation in the 
underlying theory. Studying the angular distribution can therefore be considered as a tool to discriminate between 
different scenarios beyond the SM [H, Moreover, in case of low statistics, studying the angular distribution can 
help us to veto the background Notice, however, that with this method only the absolute values of Al and 
An can be derived and no information on the relative phase of Al and Ar can be extracted. Whereas the relative 
phase of Al and An carry valuable information on the sources of CP-violation in the underlying theory. In this 
letter, we show that if in addition to the angular distribution of the final particles in the LFV muon decay, we also 
measure their polarization, we will be able to extract the phase of A* l Ar. Remembering the fact that the state-of- 
the-art LHC experiment will most likely not be able to measure these phases and for measuring such phases a more 
advanced collider, ILC, is proposed the possibility of measuring these phases by muon decay experiments seems 
more exciting. In section II, we show that by studying the angular distribution of transversely polarized positrons 
and photons we can extract the relative phase of Al and An- In section III, we discuss the possibility of extracting 
the same information by studying the angular distribution of the final positrons produced in n + — ► e + e~e + . We then 
compare the two methods and discuss the advantages and disadvantages of each one. We summarize our results in 
section IV. 



II. LEPTON FLAVOR VIOLATING RARE DECAY fj, 



Consider an anti-muon at rest [i.e., = (m M , 0, 0, 0)] which decays into a positron and a photon with definite 
spins of s e and s 7 , respectively. Using the effective Lagrangian |T]). we can calculate the fi — » &y decay rate: 



dT[/i+(P„+) - e+(P e+ ,s e+ ) 7 (P 7 ,s 7 )] m 



3 r 



dcos( 



m 

8vr 



|« + | 2 |,4l| 2 (1 + P m cos( 



) sin 



(3) 



\a-\ 2 \A B \ 2 {l-f^ cos 0) cos 2 y + ¥ fi Rc[a + a*_A* L A R e l4 ' s } sin sin f 

where P M is the polarization of the anti-muon, is the angle between the directions of the spin of the anti-muon 
and the momentum of the positron, and 9 S is the angle between the spin of the positron and its momentum. In 
the above formula, <f> s is the azimuthal angle that the spin of the final positron makes with the plane of spin of the 
muon and the momentum of the positron (to be specific to measure <j> s , the coordinate system has been defined as 
follows: z = Pe+/\p7+\ an d y — s e + x z/\s e + x z\). Finally, a+ and a_ give the polarization of the final photon: 
e» = (0,a+ - a_)i, 0)/a/2 with y/\a+\ 2 + |«- 1 2 = 1. 

Summing over the spins of the final particles, we arrive at the well-known formula shown in Eq. @ which does not 
contain any information on the relative phase of Al and An- Moreover, from Eq. ([3]) it is clear that in order to be 
sensitive to the phase of AlA* r the combination a + a*_ sin0 s should be nonzero. Remember that a_ = and ct + = 
respectively correspond to positive and negative helicities. On the other hand, sin0 s = corresponds to either a 
right-handed positron (for S = 0) or to a left-handed positron (for S = 7r). Thus, in order to extract the relative 
phase of Al and Ar from fi + — ► e + 7 we have to study the final positrons and photons whose spins are not parallel 
to their momenta. 

Let us now consider the CP conjugate of the same process. It is straightforward to prove that the partial decay 
rate of the CP conjugate process, dT/dcos9, is given by replacing Al — > A* L and An — > A* R . In other words, we 
obtain 

~t~~~~^ ~ t^TJ = ^P M Im[a + a*e*]Im[^^]sm0sin0 s . (4) 
dcos0 dcos0 47r ^ 

As expected the difference is given by the imaginary part of AlA* r . Eq. ^ shows that if we can run the experiment 
both in the muon and anti-muon modes, we will be able to derive Im[A^A|j] even without studying the angular 
distribution of the final lepton: 

-dcos0- / -dcos0 = — ^P„Im[a+ale^»lsin0 s Im[A L Ay. 

dcos0 J d cos 9 8 
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At first sight, it may seem that the above relation is at odds with the generalized optical theorem Q which states that 
the total decay rate of a particle and an anti-particle should be equal. Notice, however that we have not summed over 
the final spins so the integrals on the left-hand side do not give the total rate of ^ — ► ej. In fact, the above equation 
shows that summing over the spin of the photon and/or the positron the difference vanishes, as expected from the 
generalized optical theorem. The effect is maximal for linearly polarized photons (i.e., a_ = ±a+ = l/\/2) and for 
the final leptons polarized in the direction perpendicular to the direction of the spin of muon and the momentum of 
the final lepton (i.e., 9 S — n/2, cf> s = tt/2). 

If we have only the anti-muon mode available (or only the muon mode available), we can still extract IulL4lAJj] by 
studying the angular distribution of the final leptons. Notice that 



- 1 ' 2 dV 



■dcos 6 



l ' 2 dT 



1/2 



d cos i 



■dcos 6 



1/2 



- -dcos6> = — ^1 

a cos 6 8tt 



V3 N 



Re[a + a*_e l ^AlA R }sm6 s , (5) 



which shows that sensitivity to Im [.A£.Ar] is maximal again for linearly polarized photons (i.e., a_ = ±a 
and leptons polarized in the direction perpendicular to the direction of the spin of muon and the momentum of the 
final lepton (i.e., 9 S = tt/2,4> s = tt/2). 

Measuring the transverse polarization of the final lepton is feasible. In fact, this technique has long been employed 
to measure the Michel parameters ML Measuring the linear polarization of photon at energies of ~ 50 MeV also seems 
practical Q • As recently shown in jiol ] equipping the experiments with photon polarimeters can have implications for 
studying the radiative muon decay, too. 



III. THREE BODY DECAY ^+ 





FIG. 1: Penguin diagrams contributing to /i 4 
terms in Eq.[T] 



The vertices marked with boxes are the LFV vertices from interaction 



The effective Lagrangian in Eq. [T]can give rise to /i + — * e + e~e + through the penguin diagrams shown in Fig. [T] 
The penguin diagrams are not the only diagrams that contribute to the decay mode /i + — > e + e~e + (see the appendix 
for more details). However, we expect the contribution from the penguin diagrams to be dominant because for the 
final lepton with energy close as we will see momentarily, the photon in the penguin diagram can go almost on 

shell, resulting in an enhancement by hfm^/m,). 



In fact, we expect about 90 % of the three-body /z + 



decays to result in a lepton with energy ~ m^/2. 



In this section, we show that if we measure the spin of the final lepton with energy as well as the angular 

distributions of the final particles, we can extract information on the relative phase of Al and Ar. 

Consider an anti-muon at rest with a spin at the (x, z) plane which makes an angle of 9 with the z-axis: 



Pn+ = 0,0,0) v + 



/m u (— sin - , cos - , sin - , — cos - ) 



(6) 



where and are respectively the four-momentum and the Dirac spinor of the anti-muon. Suppose the anti-muon 
decays into an electron and two positrons with the following momenta: 



P ef =(E l ,0,0,JEl 



P e + = (1?2, |.P e +|sin0 e cos(/>, |P + |sin0 e sin<£, |P + |cos0 e 



(7) 
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FIG. 2: This figure schematically depicts the direction of the momenta of the final particles in the LFV decay fi + — > efe e\ 
relative to the spin of the anti-muon in its rest frame. 



where E 2 = |-P e +| 2 + m 2 . The above angles are illustrated in Fig. [2j The four-momentum of the final electron is 

determined by the energy-momentum conservation. The diagrams contributing to [i + — > e^e^ej are shown in Fig.[T] 
The amplitude corresponding to diagram fT]-a) can be written as 

Uel v v e+ l -^v^a^{A L P R + A R P L )v e+ q a , (8) 



where q a is the four-momentum of the virtual photon in the penguin diagram: q = P u + — P + . Combining Eqs. (pi) 

and ([7]), we find q 2 = m 2 + m 2 — 2Eim^. As a result, in the limit E\ — > q 2 <C m 2 and the amplitude in 

Eq. {5} is considerably enhanced. That is while the propagator of the virtual photon appearing in Fig. (jTJ-b) is given by 
l/(P fl + —P e +) 2 ~ so, in the limit E\ — > the effect of diagram (JT]-b) can be neglected in comparison to that 

of diagram (JT]-a). Moreover, in this limit, the effects of the LFV terms other than the term in Eq. ([T]) are lower at least 
by a factor of m^/m 2 and can be also neglected (see the appendix for more details). Let us define dT Max jd cos 6 'd<f> 
as partial decay rate of fi + into a positron with energy close to m M /2 and spinor v e + = \f2E\(Q, d e , c e , 0) T : 

dT Max rE ma „ , m „/2 ^ 



dcosOdcj) s f-^ s J m ^/2~AE dE 2 dE 1 dcos 

where AE <C wi M and E max ~ m^/2 — 4m 2 /m^. Notice that we have integrated and summed over the energies and 
spins of the pair of and e~ but not over those of e\ . It is straightforward to show that 

- " [l^rf |c e | 2 (l+P^cos60 + |A«| 2 K| 2 (l-P^co S ^ (10) 



d cos 6d<j) 768tt 4 



+P M sinfl (cos(20)Re[A R A* L d e c* e } + sm(2<f>)lm[A R A* L d e c* e })] log 



4m 



2 



i^AE ' 



where, as shown in Fig. [2j ^ is the angle between the spin of the anti-muon and the momentum of e\ and <f> is the 
azimuthal angle of the momentum of e J measured from the plane made by the spin of /i + and the momentum of e+ 
[see Eq. ©]. 

After integrating over and cos# and summing over the spin of e\ (i.e., summing over states c e = 1, d e — and 
d e = 1, c e = 0), we will arrive at the familiar formula in the literature (e.g., see |ll|). However, in this case the 
information on the phase of A R A* L will be lost. In order to extract this phase, we have to be able to measure the spin 
of e\ as well as the direction of the momenta of the final states relative to the spin of the anti-muon. Let us now 
define the following ratio 



J_ 1 dcosfl[/ dcosed0 sgn(tan ( 
tfdcoseVZ" i &sgn(tan(0 + tt/4))#] 



~~ TTH ~< r 27r dT Max ~ 7777777' 
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Notice that sgn(tan^) in the integral is equal to ±1 depending on the quadrant that <j> belongs to. In principle, if the 
polarization of the anti-muon is large (i.e., P M is about 100%), this ratio can be measured in the lab. Using Eq. (jTT 
we can show that 

= Jm[A R A* L d e 



Re[A R A* L d, 



which directly gives the phase of A R A* L for a transversely polarized positron, d e — c e = l/y/2. 
Now let us compare the advantages and disadvantages of each decay mode. In general, we expect 

yr +e + V - fp°s(4) - t] - °-°° 61 - 

Br(^+ -» e+7) 3n 1 ml' 4 J 

Thus, measurement of /i + — > e + e~e + will suffer from a higher statistical uncertainty. On the other hand, to extract 
the relative phase of Al and A R by studying ji + — > e + j in addition to measuring the spin of e + , it is necessary to 
measure the spin of 7, too. Whereas in the case of [i + — ► e + e~e + , one has to measure only the spin of the final 
positron with energy close to m M /2. 

As is well-known in the case of a three-body decay mode, we can have CP- and T-odd observable quantities, even if 
the spins of the final particles are averaged over. However, the above discussion shows that if the effective Lagrangian 
([1]) is the only source of LFV, once we average over the final spins, the CP- and T-odd effects will disappear. In fact, as 
shown in [1, [T2|, if the four-fermion LFV terms listed in the appendix are also present, the CP- and T-odd effects will 
persist even after averaging over the final spins. However, we generally expect these effects to be suppressed roughly 
by a factor of Cim fl /[AL jR log(4mg/m IJ ,AE)] compared to the effect we have discussed in the present paper. Notice 
that the two effects are sensitive to different combinations of the CP-violating phases and can be thus considered as 
complementary. 



IV. CONCLUDING REMARKS 



In this paper, we have suggested a new method to derive information on the sources of CP-violation beyond those 
in the SM. The method is based on studying the rare LFV decay of polarized muons. We have performed our analysis 
within a general effective LFV Lagrangian so our results apply to any beyond SM scenario that violates the lepton 
flavor by adding new particles at energies higher than the electroweak symmetry breaking scale (e.g., supersymmetry, 
large extra dimensions). 

We have first studied the LFV rare decay fj, — > cy and shown that provided that the polarization of the final 
particles are not parallel to their momentum (e.g., if e and 7 are respectively transversely and linearly polarized), by 
studying the angular distribution of e and the photon relative to the polarization of [i, we can extract information 
on the CP-violating phases. We have also shown that if both muon and anti-muon modes are available, the same 
information can be derived by comparing T(/i + — > e^7„) and T(fi~ — > e~7„) where the subscript tt indicates that the 
spin of the particle is not parallel to its momentum. 

We have also studied the fi + — > e~ej decay where is defined to be the more energetic positron. The amplitude 
of /i+ — > e\e~e^ is severely enhanced if the energy of e\ is close to m M /2. Thus, we expect the majority of e\ to have 
energies close to ra^/2. We have focused on decays with such kinematics and proposed a new method for extracting 
information on the CP-violating phases which is based on studying the angular distribution of the final particles. We 
have shown that with transversely polarized e\ one can extract information on a combination of the CP-violating 
phases that is impossible to achieve if e\ with helicity ±1 is employed or if the final spins are averaged over. Notice 
that in this method measuring the spin of only one of the final particles (i.e., et) will be enough. We have discussed 
the differences and synergies between this method and the one discussed in [3, [jjj • 



Appendix: LFV effective Lagrangian 



In the appendix, we discuss possible LFV operators that appear by integrating out the heavy states within theories 
such as the MSSM and show that the effect of Eq. |T]) on rare LFV muon decays is dominant. 

In the literature (see e.g., [Hj|), ^ has been shown that in the context of the MSSM, integrating out the heavy 
supersymmetric states the effective LFV Lagrangian of the electron- muon system will, in addition to Eq. ([T]), contain 



e a faq 2 la (B L P L + B R P R )e + H.c, 



(12) 
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where e a is the photon field, q is the momentum of the photon, Pl (Pr) is left (right) projection matrix and Bl and 
Br are couplings with dimension of [mass] -2 . This effective term will have no impact on \i — > simply because 
for on-shell photon (q 2 = 0), this term vanishes. However, in general it can contribute to /i + — > e + e~e + through a 
penguin diagram. Notice that unlike the case of Eq. ([5]), in this case the penguin diagram does not diverge as the 
photon propagator goes on-shell. Moreover, for most of the parameter space of the MSSM B^ jtm^ <C A^ r so the 
effect of Eq. (fT2|) is further suppressed. As a result, for calculating Y Max [defined in Eq. |9])] we can neglect the effect 
ofEq. CED). 

The effective LFV effective Lagrangian will also contain the following four-fermion terms that can in principle 
contribute to /i + — > e + e~e + : 

H = C\(p, R e L ){e R e L ) + C 2 {p<LeR){e L e R ) 

+ C 3 (p, R 'y^e R )(e R 'y fl e R ) + C4(p,L^e L )(e L ^e L ) 

+ C h {jl R ^^e R ){e L ^^e L ) + C 6 {fi L ^e L ){e R ^^e R ) + YL.C. ( 13 ) 

Again, we expect that the effect of the above four-fermion terms on Y max [see Eq. ([9])] to be negligible compared to 
terms in Eq. ([1]). That is because, unlike the penguin diagrams in Fig. [H the diagrams corresponding to the above 
four-fermion interaction terms do not diverge for E\ — > Moreover, for most of the parameter space of the 

MSSM the four-fermion effective couplings are small; i.e., Ciiri^ <C Al,r. 

One should notice that the most general effective LFV Lagrangian, in addition to the terms discussed above, 
contains extra terms. For example, it is possible to have terms such as 

e a ^fjip a lp{D L P L + D R P R )eF^ + H.c, 

where p a is the four-momentum of the electron. However, studying the effective LFV Lagrangian in its most general 
form is beyond the scope of this letter. 
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